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1 Introduction 

We know already the following properties of LH media: 

1. €r and /ir are frequency dispersive. 

2. er and [ir are negative over a similar frequency band. 

3. The tryad {E, H, k) is left-handed. 

4. The index of refraction is negative. 

From the past lectures, we know that these materials can be realized by a succession of wires 
and rods: 

• Periodic arrangement of rods: realizes a plasma medium with negative 6^ over a certain 
frequency band. The model for the permittivity is: 



l^r = l--, , (2) 



• Periodic arrangement of rings (split-rings) realizes a resonant jir modeled as 

where F is the fractional area of the unit cell occupied by the interior of the split-ring 
{F < 1). 

In the lossless case (7e = 7m = 0), we can rewrite these two relations as: 

>..= "'-r°-f"' -(1-^)4^. (3b) 



where cj^ = cjq/VI — F > ouq. Therefore: 

k =u; efi= -{1- F) f 2 • W 

Upon identifying the regions where and fir change signs, we can immediately get the 
relation for k: 
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Section 2. Argument on n < 
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The region u G [c<;o,c^6], which also corresponds to e^- < and jir < 0, corresponds to 
positive, which means k real. Therefore, there is propagation in this band, but not in the 
adjacent ones. 

It may still not be clear that k is negative, even if we write 

k = uj^^/eji = u'^^eliil ^Jer^ir = k^n . (5) 
A demonstration of the fact that n is negative follows. 

2 Argument on n < 

2.1 Complex Poynting theorem 

We shall first recall the derivation of the complex Poynting theorem and the signification of the 
various terms. 

We start from Maxwell's curl equation 

V X iuB , (6a) 

V xH = -iuoD + J . (6b) 

Upon multiplying Eq. (6a) by and substracting the complex conjugate of Eq. (6b) mul- 
tuplied by E we get: 

. V X E-E V xH"" = V '{ExH'') 

= iujB . H"" - iuoD'' - E - J"" - E 

= i(jj[B ' S'' - E ' D""]- E ' J"" . (7) 

Upon rewriting, we get: 

. = V • (^ X ^^) + iu;[E D'' -B'H'']. (8) 

On the right-hand side of the equation, the first terms corresponds to the divergence of 
Poyting power, which is therefore positive. The second term relates to the complex EM energy, 
and is therefore also positive. Consequently, the left-hand side term must also be positive, and 
actually corresponds to the power supplied by J to the volume. 

We shall use this result hereafter. 
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2.2 ID wave equation 

For the sake of simplification, let us work with a ID problem. The wave equation 

V^E{f) + k'^Eif) = -iu;fiJ{f) , (9) 

is rewritten with 

E{f) =zE{x), (10a) 
J(f) =zjo6{x - xo) , (10b) 



to yield 



q2 

^E{x) + k'^Eix) = -iufijoSix - xo) . (11) 



dx'^ 

The solution to this equation is 

^(x) = c^e^^l™l, (12) 
where a needs to be determined. From Eq. (12), we write: 

1. First derivative: 

^^ = aik-^\x-xo\ e^'^l^-^ol. (13) 

2. Second derivative: 

=aik ^k-^ol e^^l^-^ol +c.(-A:2)(^|x-xo|)^e^^l---ol 

= -ak'^ e^^l^-^ol + 2iakS{x - xo) . (14) 



Therefore: 



d'E{x) 



dx'^ 

Comparing Eq. (11) to Eq. (15), we get 



+ k E{x) =2iak6{x — xq) = 2ia ko n5{x — xq) . (15) 



= = ^ 5 (16) 

so that finally the solution is: 

E(x) = ^ e'^l---°l . (17) 

2 Tt 

If we now compute the power supplied by the current J to the volume: 

P = -i / E-rdV = ^^ >0. (18) 
2 Jv 4 n ^ ^ 

The source must, on average, do positive work on the field. Yet, in LH regime, we have 

fir < so that we must have n < as well. 
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Section 3. Dispersion relations 



Finally, we can also write the E field as: 

E{x,t) = joe^^^°'''''"''°'"^'^ . (19) 

Thus, plane waves appear to propagate from — oc and +oc to the source, seemingly running 
backward in time. Yet, the work done on the field is positive so clearly the energy propagates 
outward from the source. 



3 Dispersion relations 

At this point, we know that n < and k < 0. The difference between phase and group velocity 
can be directly seen on the dispersion relation diagram. 

v^=^ , (20a) 




(20b) 



• Free-space: k = ojy/ejl where e = cte and fi = cte. 

• Metamaterial: 



Let us take the following models 



er =1 — 



Mr =1 



(21a) 
(21b) 



3.1 Lossless case {je = 7m = 0), = 

We rewrite 
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(22b) 
(22c) 



and plot the relations with 
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k surface (Gamma= 0) 



3.1 Lossless case {je = 7m = 0)7 



Other cases follow. 
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3.1 Lossless case {je = 7m = 0)7 
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Plotting all the 3D curves on the same scale: 




(a) 7 = rad/s. (b) 7 = 10e6 rad/s. 



k surface (Gamma= 1 00000000) , ,0'° k surface (Gamma= 1 000000000) 




(c) 7 = 10e7 rad/s. (d) 7 = lOeS rad/s. 



k surface (Gamma= 1.000000e+10) ,10'° k surface (Gamma= 1.000000e+11) 




(e) 7 = 10e9 rad/s. (f) 7 = lOelO rad/s. 

Figure 1: Dispersion relation (real k) for various values of losses 7. The models are: 
Cr = 1 - (joI / {uj'^ + i^uj). fir = ^ - {^mp " ^mo) I {^'^ " ^mo + ^7^)? wliere CJp = 30e9 
rad/s, cjmp = 20e9 rad/s, cjmo = 5e9 rad/s. 



12 3.1 Lossless case (7e = 7m = 0), 



For simplicity, we can study the lossy case for 7e = 7m amd uUmo = (although we don't 
really simulate the same medium, the fundamental behavior is similar, and simpler to carry out 
mathematically) . 

The model therefore reads: 



er = fir = o-T^^ . (23) 



We compute: 



u'^ — (jjp + iju 



v2, ,2 



(24a) 



The real part is given by: 

Losses have the effect to lower the plasma frequency to 

^^ = \l^l-l^- (26) 

In addition, we also see that if 7 is very large, the plasma effect will completey dissapear 
(cf. dispersion relation for 7 = lOelO rad/s). 



